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Chiral Symmetry in Strongly Interacting Matter 

From Nuclear Matter to Phases of QCD 
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Physik- Department, Technische Universitdt Miinchen 
D-85747 Garching, Germany 

This is a brief summary of topics that were presented as lectures within the programme 
"New Frontiers in QCD 2010" at the Yukawa Institute of Theoretical Physics in Kyoto. The 
basic subject is phases and symmetry breaking patterns as they emerge from the approximate 
chiral symmetry of QCD. Part I focuses on the QCD interface with nuclear physics via chiral 
effective field theory. This includes nuclear thermodynamics and, in particular, constraints 
for compressed and hot baryonic matter provided by the density and temperature dependence 
of the chiral condensate. Part II explores aspects of the QCD phase diagram using a non- 
local improved version of the Polyakov - Nambu - Jona-Lasinio (PNJL) model. A prominent 
O j feature of such an approach is the occurence of a dynamical entanglement between chiral and 

deconfinement crossover transitions. Comparisons with available results from lattice QCD 
thermodynamics will be made. 

^ §1. Part I: Nuclear Chiral Dynamics 

c| Chiral effective field theory is the low-energy realization of QCD in the meson 

and single-baryon sectors. It is also an appropriate framework for dealing with the 
nuclear many-body problem in terms of in-medium chiral perturbation theory0 In 
this approach, chiral one- and two-pion exchange processes in the nuclear medium 
are treated explicitly while unresolved short-distance dynamics are encoded in con- 
tact interactions. Three-body forces emerge naturally and play a significant role 
in this framework. The pion mass rn^ , the nuclear Fermi momentum pp and the 
mass splitting Ma — Mjy — Zm^ between the Z\(1232) and the nucleon are all com- 
parable "small" scales that figure as expansion parameters. The relevant, active 
degrees of freedom at low energy are therefore pions, nucleons and A isobars. Two- 
pion exchange interactions produce intermediate-range Van der Waals - like forces 
involving the large spin-isospin polarizablity of the individual nucleons. The Pauli 
principle plays an important role acting on intermediate nucleons as they propagate 
in two-pion exchange processes within the nuclear medium. 
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1.1. In-medium chiral perturbation theory 

A key ingredient of this approach is the in-medium nucleon propagator, 

-M N + ie)- 1 - 2tt(7^ + M N ) S(p 2 - M 2 N ) 9(p°) 6{p F - \p |) , (1-1) 

that takes into account effects of the filled nuclear Fermi sea (with Fermi momentum 
Pf). These propagators enter in the loop diagrams generating the free energy density. 
Thermodynamics is introduced using the Matsubara formalism. The free energy 
density is then computed as a function of temperature T at given Fermi momentum 
pp. The loop expansion of this in-medium chiral perturbation theory is connected 
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in a one-to-one correspondence to a systematic expansion of the energy density in 
powers of p^, with expansion coefficients expressed as functions of the dimensionless 
ratio m^/pF. 

Calculations have been performed up to and including three-loop order in the 
energy density. A very limited set of constants associated with NN contact terms 
is fixed by reproducing e.g. the binding energy per nucleon in equilibrium nuclear 
matter at T — 0. Then this scenario leads to a realistic nuclear matter equation 
of state 1 ) with a liquid-gas first order phase transition and a critical temperature 
of about 15 MeV (see Fig.l), close to the range of empirical values extracted for 
this quantity. Incidentally, this nuclear liquid-gas transition is so far the only well 
established part of the phase diagram for strongly interactingmatter at finite density, 
with the possible exception of the chemical freezeout curved deduced from hadron 
production in heavy-ion collisions. 

The three-loop truncation in the chiral expansion of the pressure in powers of 
the Fermi momentum, pp, implies that these calculations can be safely trusted up 
to about twice the density of normal nuclear matter, i.e. for pp < 0.3 GeV << 4^/^ 
with f n = 0.09 GeV the pion decay constant in vacuum. 




Fig. 1. Equation of state (left) and phase diagram (right) of symmetric (N = Z) nuclear matter 
from in-medium chiral effective field theory0'Hl The critical point for the first-order liquid-gas 
phase transition appears at T c ~ 15 MeV. The equlibrium point at T = and baryon density 
po ~ 0.16 fm -3 is fixed by a single parameter, the strength of an NN contact interaction. All 
remaining input is pre-determined by known pion-nucleon interactions in vacuum. 

Figure 2 shows (T, p) phase diagrams for isospin-asymmetric nuclear matter, 
outlining the evolution of the critical point and the gas-liquid coexistence region 
for increasing proton fractions Z/A. The first order phase transition stops exist- 
ing at Z/A — 0.05, at which point the liquid component disappears and neutron 
matter begins to be realized as an interacting Fermi gas. It is important to note 
that this behaviour is almost entirely driven by the isospin dependence of the in- 
medium two-pion exchange interactions. Typical ingredients of one-boson exchange 
phenomenology (such as rho and sigma bosons) do not appear in this approach and 
are replaced by explicit 2tv exchange mechanisms. 

Extensions of this framework to finite nuclei, using the energy density func- 
tional formalism, have been applied successfully 4 * to calculate ground state proper- 
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Fig. 2. Phase diagram for isospin-asymmetric nuclear matter from in-medium chiral effective field 
theory 2 ' plotted in the (T, p) plane with p — p p + p n as function of the proton fraction Z/A. 
The trajectory of the critical point is also shown. 

ties (binding energies, mean square radii, deformations and systematics along isotopic 
chains) throughout the nuclear chart, from ie O to Pb isotopes. Other examples of 
detailed nuclear structure phenomena such as the anomalously weak 14 C beta decay 
transition, with its lifetime of almost six thousand years that enables radiocarbon 
dating, can be understood in terms of the density dependent chiral effective interac- 
tion including a pronounced three-body contribution. 5 ) 

1.2. Density and temperature dependence of the chiral condensate 

In a nuclear equation of state based on chiral dynamics the pion mass enters 
explicitly (or, equivalently, the quark mass m q according to the Gell-Mann - Oakes 
- Renner relation, rn^f^ — —m q (tpt(j)). Equiped with such an equation of state one 
can now ask the following questions: how does the chiral (quark) condensate (V^) 
extrapolate to baryon densities p exceeding those of normal nuclear matter, and 
how does it evolve as a function of temperature T in a nuclear environment? The 
answer, first at T — 0, is based on the Hellmann-Feynman theorem applied to the 
Hamiltonian density of QCD in the nuclear ground state \Jf), 



where (ipip)o is the chiral condensate in vacuum. The energy density £ (normalized 
so that it vanishes at p 0) is expressed at given baryon density p as a function of 
the pion mass. In-medium chiral effective field theory gives the following result: 



9Uqcd 
dm q 



\V) = {W)o 1- 



1 dS(m^p) 
SI dml 



(1-2) 




+ ... - 



P d 
ft dm* 




The second term on the r.h.s., with its leading linear dependence on density, is 
the contribution from a free Fermi gas of nucleons. It involves the pion-nucleon 
sigma term cr/v — 0.05 GeV and additional non-static corrections. The third term 



4 



W. Weise 



includes the pion mass dependence of the interaction energy per nucleon, Ei nt /A. 
This term features prominently the two-pion exchange interaction in the nuclear 
medium including Pauli principle corrections, and also three- nucleon forces based on 
two-pion exchange. 
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Fig. 3. Left: density dependence of the chiral condensate at temperature T — in symmetric 
nuclear matter P Dot-dashed curve: leading order term using a at = 50 MeV. Upper curve: 
full in-medium chiral effective field theory result at three-loop order. Lower curve: chiral limit 
with vanishing pion mass. Right: temperature dependence of the chiral condensate starting 
from the upper curve of the left figure at T = and using finite-temperature in-medium chiral 
perturbation theoryEf 



The dashed-dotted curve in Fig. 3 (left) shows the pronounced leading linear 
reduction in the magnitude of the chiral condensate with increasing density. This 
holds in the absence of correlations between the nucleons. Up to about the density 
of normal nuclear matter, this term dominates, whereas the interaction part tends 
to delay the tendency towards chiral restoration when the baryon density is further 
increased. This stabilizing effect arises from the combination of Pauli blocking in 
two-pion exchange processes and three-body correlations, and their explicit depen- 
dence on the pion mass. This behaviour is highly sensitive to the actual value of the 
pion mass. In the chiral limit, — >► 0, with stronger attraction in the NN force at 
intermediate ranges, the trend is reversed and the more rapidly dropping condensate 
would now lead to the restoration of chiral symmetry at relatively low density. In 
essence, nuclear physics would look radically different if pions were exactly massless 
Nambu-Goldstone bosons. The drastic influence of explicit chiral symmetry break- 
ing in QCD through the small but non-zero u and d quark masses on qualitative 
properties of nuclear matter is quite remarkable. 

The additional temperature dependence of the in-medium chiral condensate is 
found by repeating the calculation, Eq.( |l-2[ ), using the free energy density, T{m^\ p, T): 

(W)p,t = 1 _ 1 dF{m«-p,T)) 

<V^)o ft dml 1 ' j 

The result shown in Fig. 3 (right) demonstrates how the chiral condensate "melts" 
with increasing temperature and approaches the linear density dependence of the 
non-interacting Fermi gas at a temperature of about 100 MeV. The effect of thermal 
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pions on the temperature dependence of the quark condensate must still be added; 
it amountdlMIl to lowering the magnitude of the condensate by another 5% of its 
vacuum value at T ~ 100 MeV. However, no tendency towards a first order chiral 
phase transition is visible, at least up to baryon densities about twice that of normal 
nuclear matter and up to temperatures T ~ 100 MeV. In that range it appears 
that the relevant fermionic quasiparticle degrees of freedom are indeed nucleons 
rather than quarks. This should be kept in mind as an important constraint for the 
discussion of the QCD phase diagram at finite baryon chemical potential. 

§2. Part II: Phases of QCD - chiral and deconfinement transitions 

Confinement and spontaneous chiral symmetry breaking in QCD are governed 
by two basic symmetry principles: 

i) The symmetry associated with the center Z{3) of the local SU(3) C color gauge 
group is exact in the limit of pure gauge QCD, realized for infinitely heavy quarks. 
In the high-temperature, deconfinement phase of QCD this Z{3) symmetry is spon- 
taneously broken, with the Polyakov loop acting as the order parameter. 

ii) Chiral SU(Nf)ft x SU(Nf)L symmetry is an exact global symmetry of QCD 
with Nf massless quark flavors. In the low-temperature (hadronic) phase this sym- 
metry is spontaneously broken down to the flavor group SU (Nf)y (the isospin group 
for Nf = 2 and the "eightfold way" for Nf = 3). As a consequence there exist Nj — 1 
pseudoscalar Nambu-Goldstone bosons and the QCD vacuum hosts a strong quark 
condensate. 

There is no principal reason why the chiral and deconfinement transitions should 
be closely connected. Nonetheless this appears to be the case in lattice QCD com- 
putations with almost physical quark masses. Examples of results from lattice QCD 
thermodynamics^ with 2+1 flavors (at zero bar yon chemical potential) are shown 
in Fig. 4. More recently, improved computationsI2rlIII (see Fig. 5) indicate a shift 
of the chiral transition to somewhat lower temperatures, consistent with earlier lat- 
tice simulationsl!2f Nevertheless, given the smoothness of the crossover, the chiral 
and deconfinement transitions appear to be linked and overlapping at zero chemical 
potential, \i = 0. 

2.1. Chiral condensate and Polyakov loop 

The order parameter of spontaneously broken chiral symmetry is the quark con- 
densate, (V^)- The disappearence of the condensate, by its melting above a transi- 
tion temperature range T c , signals the restoration of chiral symmetry in its Wigner- 
Weyl realization. The transition from confinement to deconfinement in QCD is like- 
wise controlled by an order parameter, the Polyakov loop. A non- vanishing Polyakov 
loop ^ reflects the spontaneously broken Z(3) symmetry characteristic of the decon- 
finement phase. The Polyakov loop vanishes in the low-temperature, confinement 
sector of QCD. 

One observes that the chiral and deconfinement transitions as shown in Fig. 4 
are not phase transitions but smooth crossovers, so there is no critical temperature 
in the strict sense. It is nevertheless useful to define a transition temperature band 
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Fig. 4. Lattice QCD results (Nf = 2 + 1) of chiral condensate (left) and Polyakov loop (right) as 
functions of temperature 1^ The grey bands mark the crossover regions and give an impresion 
of the transition temperature range. Different data sets correspond to different number N T of 
lattice points along the Euclidean time axis. 
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Fig. 5. Stepwise improved lattice simulations of the chiral condensate from the Hot QCD collabo- 
ration^ in comparison with the Wuppertal-Budapest results!^ EH 



around the maximum slope of either the condensate (ipip)(T) or the Polyakov loop 
<Z>(T). 

Two limiting cases are of interest in this context. In pure gauge QCD, corre- 
sponding to infinitely heavy quarks, the deconfmement transition is established as 
a first order phase transition with a critical temperature of about 270 MeV. In the 
limit of massless u and d quarks, on the other hand, the isolated chiral transition ap- 
pears as a second order phase transition at a significantly lower critical temperature. 
This statement is based on calculations using Nambu - Jona-Lasinio (NJL) type 
models which incorporate the correct spontaneous chiral symmetry breaking mech- 
anism but ignore confinement. The step from first or second order phase transitions 
to crossovers is understood as a consequence of explicit symmetry breaking. The 
Z(3) symmetry is explicitly broken by the mere presence of quarks with non-infinite 
masses. Chiral symmetry is explicitly broken by non-zero quark masses. But the 
challenging question remains whether and how the chiral and deconfmement tran- 
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sitions get dynamically entangled in just such a way that they finally occur within 
overlapping transition temperature intervals. 

2.2. The non-local PNJL model 

Insights concerning this issue can be gained from a model based on a minimal 
synthesis of the NJL-type spontaneous chiral symmetry breaking mechanism and 
confinement implemented through Polyakov loop dynamics. This PNJL 
is specified by the following action: 

w3=l/T r i V 

S = l dr J v d 3 x[^d T 4>-U(4>,?p\<j))\ --U{$,T) . (2-1) 

It introduces the Polyakov loop, 

$ = AT 1 Tr exp(i0/T) , (2-2) 

with a homogeneous temporal gauge field, <j) — 03 A3 + <j)%\% E SU(3), coupled to the 
quarks. The dynamics of ^ is controlled by a Z{3) symmetric effective potential ZY, 
designed such that it reproduces the equation of state of pure gauge lattice QCD 
with its first order phase transition at a critical temperature of 270 MeV. The field 
(f) acts as a potential on the quarks represented by the flavor doublet (for Nf = 2) or 
triplet (for Nf — 3) fermion field tp. The Hamiltonian density in the quark sector is 

U = -i^{a • V + 74 m - (j))^ + V(V>, ^) , (2-3) 

with the quark mass matrix rh and a chiral SU{Nf)i, x SU(Nf)n symmetric inter- 
action V. 

Earlier two-flavor versions of the PNJL modeEU-EH have still used a local four- 
point interaction of the classic NJL type, requiring a momentum space cutoff to 
regularize loops. A more recent versionjMf using a non-local interaction, does not 
need any longer an artificial NJL cutoff. It generates instead a momentum de- 
pendent dynamical quark mass, M(p), along with the non- vanishing quark conden- 
sate. A further extension to three quark flavoraiZI includes a U(1)a breaking term 
implementing the axial anomaly of QCD. This term is constructed as a non-local 
generalization of the Kobayashi-Maskawa-'t Hooft 3x3 determinant interaction. 

A basic relation derived in this non-local PNJL model is the gap equation gen- 
erating the momentum dependent quark mass self-consistently together with the 
chiral condensate. Its form (written here for the two-flavor case) is reminiscent of 
the corresponding equation emerging in Dyson-Schwinger approaches to QCD: 

M( P ) = », + «¥/|«- ! ) ? *j, (24) 

where mo = m u = is the current quark mass, G is a coupling strength of dimen- 
sion (length) 2 and C(q) is a momentum space distribution, the Fourier transform of 
which represents the range over which the non-locality of the effective interaction 
between quarks extends in (Euclidean) 4-dimensional space-time. This distribution 
is normalized as C(q = 0) = 1. Fig. 6 (left) shows C(p) in comparison with the 
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"classic" NJL cutoff by a step function in momentum space. This translates into 
the dynamical quark mass of Fig. 6 (right). The non-local approach now permits to 
establish contacts with the high-momentum limit of QCD with its well established 
behaviour M(p) oc —a s (p)(ipip)/p 2 at p oo. At p < 1 GeV the distribution C(p) is 
designed to follow lattice QCD and Dyson-Schwinger results, or it may be deduced 



from an instanton liquid model (see Ref. 17) for more details). The non-local PNJL 



model can indeed be "derived" from QCD as demonstrated in Refs. 18), 19). 
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Fig. 6. Left: distribution C(p) representing the non-local effective interaction between quarks. The 
momentum space cutoff function of the "old" local NJL model is also shown. Right: dynamical 
quark mass (solid curve) resulting from the gap equation of the non-local model. Lattice QCD 
data extrapolated to the continuum and chiral limits— j are displayed for comparison, as well as 
the result deduced from an instanton liquid model. Figures adapted from Ref. [17| . 



2.3. PNJL thermodynamics 

Once the input is fixed at zero temperature by well-known properties of the 
pseudoscalar mesons, the thermodynamics of the PNJL model can now be inves- 
tigated with focus on the symmetry breaking pattern and on the intertwining of 
chiral dynamics with that of the Polyakov loop. The primary role of the Polyakov 
loop and its coupling to the quarks is to suppress the thermal distribution functions 
of color non-singlets, i.e. quarks and diquarks, as the transition temperature T c is 
approached from above. Color singlets, on the other hand, are left to survive below 
T c . This is seen by analyzing the relevant piece of the thermodynamic potential 
Q — —(T/V) InZ, the one involving the quark quasiparticles: 



oc 



+ 



/ 
/ 



d 3 p 
(2vr) 3 
d 3 p 



In 



In 



1 + 3$ e - {E p-^ /T + 3<T e - 2 (£p-<")/ T + e -KE P -y)/T 

1 + 30* e -(E p +ii)/T + 3< p e -2(E p +»)/T + e -3(£ p +M)/T 



(2-5) 



This suppression of color non-singlets in the hadronic phase below T c should, how- 
ever, not be interpreted as dynamical confinement. The unsuppressed color singlet 
three-quark degrees of freedom are not clustered but spread homogeneously over all 
space. At this stage, nucleons are not treated properly as localized, strongly cor- 
related three-quark compounds plus sea quarks. Thus one cannot expect that such 
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Fig. 7. PNJL model calculations of chiral and deconfinement transitions EH~E3 Left: two- flavor 
PNJL results with dashed curves showing the quark condensate in the chiral limit and the 
Polyakov loop while they are still decoupled. Coupling the quarks to the Polyakov loop back- 
ground produces shifts to the corresponding solid curves, i.e. the entanglement of chiral and 
deconfinement transitions. Right: three- flavor non-local PNJL model results including explicit 
chiral symmetry breaking by non-zero quark masses m u ,d and m s . Also shown are lattice QCD 
results 91 for 2+1 quark flavors with almost physical masses. 



a model describes properly the low-temperature phase of matter at finite baryon 
chemical potentials around 1 GeV, the domain of nuclear many-body systems. 

At zero baryon chemical potential, a remarkable dynamical entanglement of the 
chiral and deconfinement transitions is nonetheless verified in the PNJL model, as 
demonstrated in Fig. 7 (left) for the two- flavor case. In the absence of the Polyakov 
loop the quark condensate (left dashed line), taken in the chiral limit, shows the 
expected 2nd order chiral phase transition, but at a temperature way below and 
far separated from the 1st order deconfinement transition controlled by the pure- 
gauge Polyakov loop effective potential U (right dashed line). Once the coupling 
of the Polyakov loop field to the quark density is turned on, the two transitions 
move together and end up at a common transition temperature around 0.2 GeV. 
The deconfinement transition becomes a crossover (with Z{3) symmetry explicitly 
broken by the coupling to the quarks), while the chiral phase transition remains 2nd 
order until non-zero quark masses, m u ^ — 4 MeV, induce a crossover transition as 
well. 

It is interesting to compare this PNJL description of the interplay between the 
chiral and deconfinement transitions, with a scenario based on strong-coupling lattice 
QCD with inclusion of Polyakov loop dynamicsHU . Although the starting points and 
frameworks are quite different in those two approaches, recent studies 21 ' arrive at 
similar results concerning the close correspondence between chiral and deconfinement 
transition temperature ranges. 

The right side of Fig. 7 shows the three- flavor PNJL resu it© together with Nf = 
2+1 lattice data. 9) A direct comparison is clearly not appropriate but the similarity 
of the crossover transition patterns is striking, given the simplicity of the model and 
the fact that these results are derived from a mean-field approach. Further important 
steps toward systematic investigations beyond mean field are being pursuedHHf An 
important step turns out to be the incorporation of pion degrees of freedom as quark- 
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Fig. 8. Effects of mesonic quark- ant iquark modes in PNJL model calculation£3'IHI of the pressure 
(left) and the chiral condensate (right) as functions of temperature. Mean field results (MF) 
are shown in comparison with those including mesonic contributions (solid curves). Also shown 
is the chiral condensate from lattice QCD computations of the Wuppertal-Budapest groupE^ 



antiquark collective modes using RPA techniques, consistently expanding beyond 
the mean field level of the non- local PNJL modeOO . Not surprisingly, thermal 
excitations of Nambu-Goldstone bosons add to the pressure below the transition 
temperature T c , as seen in Fig. 8 (left), but these modes disappear rapidly above T c 
where the quark quasiparticles and the gluonic degrees of freedom take over. 

Pions also have a pronounced on the chiral condensate as it ap- 

proaches the transition region from the hadronic phase. The chiral crossover becomes 
significantly smoother close to T c as displayed in Fig. 8 (right). Lattice QCD re- 
sult^ with physical pion masses and realistic extrapolations to the continuum limit 
show a similar tendency. 

§3. Scenarios at finite baryon density 

3.1. Modeling the phase diagram of QCD 

Undoubtedly a prime challenge in the physics of strong interactions is the explo- 
ration of the QCD phase diagram at non-zero baryon density, extending from normal 
nuclear matter all the way up to very large quark chemical potentials n at which 
color superconducting phases are expected to occur. PNJL calculations at finite fi 
give a pattern of the chiral order parameter in the (T, /a) plane showing a crossover 
at small n that ends at a critical point. From there on a first-order transition line 
extends down to a quark chemical potential fi ~ 0.3 - 0.4 GeV at T = 0. 

The typical phase diagram that emerges from the non-local three-flavor PNJL 
model 17 ' is shown in Fig. 9. The phase to the right of the 1st order transition line, but 
below the deconfmement boundary that decouples from the chiral transition beyond 
the critical point, has been named "quarkyonic't^f and is hypothetically assumed to 
exist until superconducting phases take over at large quark chemical potential. A 
recent discussion about the possible relationship between such a hypothetical phase 
and the hadron rates produced in high-energy heavy-ion collisionaHlP may not yet be 
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conclusive but is pursued with great intensity. 

One must recall at this point that calculations and extrapolations based on 
PNJL or related approaches still have limited validity at moderate and high baryon 
densities since they do not properly incorporate baryonic degrees of freedom. They 
miss the important constraints imposed by the existence of nucleons and nuclear 
matter. This becomes evident when converting the PNJL phase diagram, Fig. 9, 
into a plot that translates quark chemical potentials into baryon densities via the 
derivative of the pressure, ps = dP/dfiB, with respect to baryon chemical potential 
\±B — 3/ig. The 1st order transition line then converts into a coexistence region along 
the density axis, as illustrated in Fig. 10 (in this case deduced from a local PNJL 
model calculator©). It becomes immediately evident that this coexistence region 
encloses what is well known as nuclear physics terrain with its transition from Fermi 
liquid to Fermi gas - but now with the "wrong" quasiparticle degrees of freedom: 
PNJL quarks rather than interacting nucleons. Consequently, interpretations based 
on the schematic phase diagram of Fig. 9 along the chemical potential axis should 
be exercised with caution. 



1.0 



T 

x c 



0.5 



deconfinement 



chiral 

crossover 



hadronic 




1 st order chiral 
Polyakov loop 



"quarkyonic" ? 



0.05 0.10 



0.15 0.20 
^[GeV] 



0.25 0.30 




700 800 



900 1000 

Pb [MeV] 



1100 1200 



Fig. 9. Three- flavor non-local PNJL model 
calculation 261 of the phase diagram in 
the plane of temperature and non-strange 
quark chemical potential. At the critical 
point the chiral crossover turns into a 1st 
order transition. The deconfinement tran- 
sition band represents the crossover range 
observed in the Polyakov loop. 



Fig. 10. Isotherms of the PNJL phase diagram 
plotted in the plane of baryon density ver- 
sus baryon chemical potential \ib — S^quark- 
This calculation is based on a local version 
of the 3-flavor PNJL model. 261 The domain 
relevant to nuclear physics is indicated as "nu- 
clear terrain" . 



3.2. The quest for the critical point 

Several further important questions are being raised in this context. A principal 
one concerns the existence and location of the critical point in the phase diagram. 
Extrapolations from lattice QCD, either by Taylor expansion^!!! around fi = or by 
analytic continuation from imaginary chemical potential)!!! have so far not reached 
a consistent conclusion. A second question relates to the sensitivity of the first order 
transition line in the phase diagram with respect to the axial U(1)a anomaly in QCD. 
This issue has been addressed in Ref . 28 ) using a general Ginzburg-Landau ansatz for 
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the chiral SU(3) effective action with axial U(l) anomaly. It was pointed out that, 
depending on details of the axial U(1)a breaking interaction, a second critical point 
might appear such that the low-temperature evolution to high density is again just 
a smooth crossover, or the first order transition might disappear altogether and give 
way to a smooth crossover throughout. In a three- flavor NJL type realization of this 
approachJH* it has recently been demonstrated how the interaction between chiral 
condensate and diquark degrees of freedom, based on the genuine U(l) a breaking six- 
point vertex, can open up a smooth transition corridor along the chemical potential 
axis in which chiral and diquark condensates may coexist. 
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Fig. 11. Trajectory of the critcal point in the 
(T, ji) plane as function of K, the strength 
of the axial U(l) breaking interaction of 
the 3-flavor PNJL model 261 . Here K is 
the coupling strength that reproduces the 
j] mass in vacuum. 



Fig. 12. A sketch of a PNJL-inspired phase dia- 
gram with a smooth crossover corridor along 
the chemical potential axis. 



An impression of the explicit dependence of the critical point on the axial 
anomaly can be obtained using the (local) three-flavor PNJL model with inclusion of 
a U(1)a breaking Kobayashi-Maskawa-'t Hooft determinant interaction and varying 
the coupling strength K of this interaction© S It turns out that the location of 
the critical point in the phase diagram varies indeed strongly with K and may even 
disappear altogether below a certain value of K (see Fig. 11). It should be added, 
though, that this sensitivity to the strength of the axial anomaly appears to be less 
pronounced in the non-local three-flavor PNJL modeliZI . 



§4. Summary and Outlook 



From the variety of existing model calculations (including those using PNJL 
models) one might draw the presumably premature conclusion that critical phenom- 
ena occur already at a density scale not much higher than that of normal nuclear 
matter. However, these models are so far not capable of working with the proper de- 
grees of freedom around and below baryon chemical potentials ~ 1 GeV (correspond- 
ing to quark chemical potentials around 0.3 GeV). Approaching the corresponding 
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density scale from below, it is obvious that constraints from what we know about 
the nuclear matter equation of state must be seriously considered (see Fig. 12). This 
has been the essence of the first part of this presentation. Implementing such nuclear 
matter constraints (that are well reproduced using the framework of chiral effective 
field theory) into a consistent picture together with the non-local PNJL model (that 
has proven to be quite successful at high temperatures) remains a major challenge 
for the future. 

Finally, in the discussion of high density, low temperature constraints on the 
QCD equation of state, recent new developments concerning the mass-radius relation 
deduced from observations of neutron stars in binarie Jlll'IHI are potentially of great 
interest in reducing the set of acceptable equations of state. 
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